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ABSTRACT 



The charge distribution of periodic electron beam pulses generates patterns 
of sub-Cerenkov radiation distinctive of the distribution of charge with in a 
bunch. Mapping the radiation pattern from different charge shapes may provide 
insight into whether charge pulse shapes can be determined from observed 
radiation patterns. The radiation patterns of Gaussian, Level, and Trapezoidal 
function were mapped by computer simulation. Near 90° to the beam, the 
radiation patterns of all three charge distributions developed an envelope 
proportional to the fourier transform of the charge bunch distribution when the 
wavelength of the emitted radiation was comparable to the size of the bunch. 
For the Gaussian function, the envelope is Gaussian, for the level function it 
is a sine function. Since the envelope for the trapezoidal function in the 
product of two sine functions it is more difficult to analyze. This work may 
provide a basis for determining the charge shape of electron beam pulses from 
the sub-Cerenkov radiation based on the radiation intensity pattern. 
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I. INTRODUCTION 



A. HISTORY 

A very faint emission of bluish-white light, what was later-on called 
Cerenkov radiation, from transparent substances placed nearby a radioactive 
sources had been observed by Mnie Curie in 1910 . Although electromagnetic 
theory had been sufficiently developed at that time to describe this 
phenomenon, the first deliberate attempt to study the phenomenon was made 
by Mallet after many years, who found that the light emitted from a wide 
variety of transparent bodies placed close to radioactive source always had the 
same bluish-white light quality, and that spectrum was continuous, not 
possessing the line or band structure characteristic of fluorescence [Ref. 1: p. 
!]■ 

P.A. Cerenkov carried out . a series of experiments between 1934 and 193S 
related to the phenomenon. In 1937 Ilya Frank and Igor Tamm proposed a 
satisfactory theory of the radiation. The experimental results and theoretical 
predictions were in excellent agreement. A complete description of Cerenkov 
radiation, as well as exact mathematical treatment is given in [Ref. 1: p. 1]. 

B. WHAT IS THE CERENKOV RADIATION ? 

A charged particle or a group of charged particles moving through a 
medium at a speed greater than the (phase) velocity of light in the same 
medium will generate a continuous spectrum of electromagnetic radiation. This 
process is referred to as Cerenkov radiation. [Ref. 2: p.9] 
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C. PREVIOUS RESEARCH AT N.P.S. 

The research on Cerenkov radiation has been conducted by Professors Fred 
R. Buskirk, John R. Neighbours and others in the Naval Postgraduate school 
since the LINAC was established. Recent work has concentrated on the 
determination of electron beam bunch profile through Cerenkov radiation. 

The work was aimed to illustrate examples of the consequences of the 
finite spatial charge distribution to show that Cerenkov radiation emitted at 
low microwave frequencies may be used to characterize the properties of an 
intense relativist; beam in air. The results were that the spatial charge 
distribution of an electron pulse, along with the beam interaction length, 
determines the Cerenkov radiation distribution as a function of frequency. An 
angular distribution of the Cerenkov radiation can, in principle, measure its 
spatial charge distribution. [Ref. 3:p. 1994]. 

D. PURPOSE 

The charge distribution of a bunch influences the pattern of radiation 
produced. Detailed calculations for several simple charge distributions have 
been carried out for the Cerenkov case [Ref. 4:pp. 1994,1996]. The object of 
this work is to assess the effect of the charge distribution of a bunch on the 
sub-Cerenkov radiation from a charged particle beam. Specifically, this means 
investigation of the short beam path region at energies below the threshold. 
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II. THEORY 



A. BACKGROUND 

A fast, charged particle may cause electromagnetic radiation in a dielectric 
medium when it moves with a constant speed. Suppose an electron to be 
moving relatively slowly through a transparent medium as shown on Fig. 1.1 
a). As the particle passes through the medium, each elemental region of the 
medium along the track will in turn receive a very brief electromagnetic pulse 
and be polarized, producing a dipole. An electromagnetic pulse is generated by 
the formation and disappearance of the dipole as the bunch passes. The 
electromagnetic pulse propagates away from the dipole source at the velocity of 
light in the medium. The distance between wave fronts is compressed as the 
velocity of the electron bunch increases. This effect is similar to the doppler 
effect for sound radiated from a moving source. If the speed of the particle is 
less than the propagation velocity in the medium, the electromagnetic field 
produced is called sub-Cerenkov radiation. [Ref. 1: pp. 3 — 4] 

If the electron is moving faster than light in the medium (see Fig 1.1 b), 
the the pulse source velocity is greater than that of the propagating 
electromagnetic wave which transport the energy. The result is a wavefront 
and where the emanated wavelets bunch together in phase. This shockfront of 
electromagnetic radiation is Cerenkov radiation. [Ref. 1: pp. 4-5] 
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B. CHARACTERISTICS OF CERENKOV RADIATION 

Fig. 2.2 shows the resultant Cerenkov radiation plane wavefront BC from 
particles traveling from A to B. Wavelets from point 0 to 5 are coherent 

along wavefront BC. The angle 0 C is given by the "Cerenkov relation": 

cos 9 C = 1/3 n, (2.1) 

where n is the refractive index of the medium and 3 is the ratio of the 

velocity of the particle to the velocity of light in a vacuum [Ref. 1: pp.4-5). 
From the above, the following is seen that: 

(1) For medium of given refractive index n, there is a threshold velocity 

$mj n = (1/n), below which no radiation take place [Ref. 1: p. 5). 

(2) For ultra-relativistic particles, for which 3 i 1. there is a maximum 

angle of emission, given by 0 max = cos''(l/n). 

(3) Radiation occurs mainly in the visible and near visible region for 
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which n > 1. In the X-ray region, above atomic resonances, n is 
usually less than 1; therefore Cerenkov radiation in the X-ray region 
is not probable [Ref. 3: p. 3252, Ref. 1: p. 6]. 




Fig. 2.2 Radiation Shockfront 

(4) In a gas, 0 C has only a slight dependence on electron velocity because 
must be close to one in order to obtain the Cerenkov condition 
[Ref. 3: p. 3249]. 

From the nature of the radiation, which has been described above, it may be 
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concluded: 



(5) Azimuthal symmetry of the radiation produces a cone of radiation. 

(6) The distribution in 9 of the radiation intensity approximate a Dirac 
^-function centered about Cerenkov angle '9 C '. 

(7) Length AB in Fig.2 must be much larger than the radiated 
wavelength in order to avoid diffraction effects. For a finite radiation 
length, the diffraction effect is significant. 

C. CALCULATED POWER OF COHERENT CERENKOV RADIATION 

The total coherent power per unit solid angle, radiated at the frequency v 
by a periodic charged particle beam in traveling a finite distance L at 
constant velocity is 



where vq is the fundamental frequency of the beam generator and u is a 
harmonic of vo, q is the charge of an individual bunch, Q is a constant and R 
is the radiation function. The bunches in the beam are assumed to be 
unchanging in shape and size as the beam travels through the medium (usually 
air) at a velocity given by v = /3 Co- The velocity of light in vacuum is Co 
and the velocity of light in the medium is c = Co/n where n is the index of 
refraction. The parameters used to describe the radiation are 



W(i/,k) = R 2 , 



( 2 . 2 ) 



Q — /iCq 2 /8^7 



(2.3) 
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which is a constant, where /i is the permeability of the medium, and the 
radiation function R is given by 

R = 2xr) sin# I(u)F(k), (2.4) 



where 9 is the angle between the direction of travel of the charged particle 
beam and the direction of propagation of the emitted radiation, I(u) is the 
diffraction function, and F(k) is the form factor. [Ref. 5: pp. 2741, 2742] 

The wave vector of the emitted radiation is k (k = u/c), and F(k) is the 
dimensionless form factor. That is, if p(r) is the charge distribution of a single 
bunch, then the Fourier components of the charge are 



00 

P(k) = 

— 00 



/?(r)e lk ' r d 3 r, 



and the form factor is defined by 



(2.5) 



p(k) = qF(k). 



( 2 . 6 ) 



For a point charge, F(k) is identically one for all values of k. 

The dimensionless beam length parameter 77 which is defined as the ratio of 
the path length of the charged particle beam to the wavelength in the medium 
of the emitted radiation, appears in the radiation function explicitly, and 
implicitly as part of the diffraction parameter defined below: 

r) = L / A. (2.7) 
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The diffraction function is 



I(u) = sin u/u, (2.8) 

where the dimensionless diffraction parameter u depends upon the Cerenkov 
angle given by cos 9 C = (n/?)‘ l , as well as upon the beam length parameter 
and the radiation emission angle: 

u — it 7] [ (n/?)" 1 - cos 9 ]. (2.9) 

The parameter t; determines the periodicity of the diffraction pattern and also 
determines its magnitude. Since R depends directly upon the strength of the 
radiation is proportional to the square of the path length of the beam. [Ref. 5: 
pp. 2741,2742] 

D. FINITE-LENGTH PATH 

If the emission region has finite length, which may be realized by passing 
the electron beam through a gas cell, the radiation propagation direction is not 
confined to a sharp Cerenkov angle 0 C , but is spread over a range of emission 
angles. For the case of gas, for which the ordinary Cerenkov angle is a few 
degrees, the broadening is asymmetric about the Cerenkov cone, and the 
radiated power displays various interference lobes. This spreading effect depends 
only on the path length in the gas and does not depend on the electron bunch 
structure or periodicity of the bunches. [Ref. 3: pp. 3246,3251] 

Radiation from the polarized molecules along the path of the electron 
beam is formally equivalent to diffraction from a single slit at angles far from 
normal incidence. The shift in the radiation peak can be understood as a 
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diffraction effect arising from the linearly varying phase of the radiation along 
the beam interaction length [Ref. 3: p. 3246]. 

Considerable power (about half in the case of a gas) appears outside the 
main lobe, and the latter is the peaked at angles much larger than the 
Cerenkov angle. The Cerenkov power radiated per unit path length is increased 
by large factors approaching 2 orders of magnitude for a gas path of finite 
length compared to one of infinite length [Ref. 3: p. 3251]. 

This increase in power is associated with the finite length of the radiating 
medium only. The effect should be approximately the same for either single or 
periodic electron bunches and should occur whether the bunches are effectively 
point charges or have significant size. 

The value of the parameter u is affected by rj and the beam energy 
through p. As either of these parameters is varied, the angular position of the 
zeros of I(u) changes. The emission threshold is defined as the situation when 
the first zero of the diffraction function enters the physical range. Applying 
this criterion, the result is that for small values of rj, a Cerenkov-like 
radiation pattern is produced even when the product n P < 1. That is, the 
threshold occurs for nP < 1 at small values of rj. As rj increases (at constant 
n P < 1), the threshold is crossed and the radiation pattern becomes more 
chaotic. A full discussion of these geometrical effects is given in Ref. 5. 

E. PERIODIC ELECTRON BUNCHES 

Cerenkov radiation was considered for periodic bunches of electrons such 
as would be emitted by a typical traveling-wave electron linear accelerator 
(LINAC). Because the distribution of intensity of Cerenkov is proportional to 
the frequency, the radiation at microwave frequencies would be low unless 
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beams are intense and bunched so that coherent radiation by many electrons 
contributes [Ref. 2: p. 1531]. 

Periodic bunches of electrons produce Cerenkov radiation at harmonics of 
the bunch frequency, while a single bunch will radiate with a continuous 
frequency distribution. The effect of bunch size is reflected in the factor of the 
Fourier transform of the charge in the bunch, which causes the radiation to 
fall off at high frequencies, while at low frequencies all charges in the bunch 
radiate in phase. 

At low frequencies such that the wavelength of the emit, ted radiation is of 
the order of the bunch size, the electrons in the bunch radiate coherently. This 
leads to large enhancement factors for typical bunches consisting of 10 8 
electrons, and in fact allows the radiation to be significant at microwave 
frequencies. [Ref. 3: p. 3246] 

At higher frequencies or harmonics, destructive interference, described by 
the Fourier transform of the charge density, decreases intensities with increasing 
frequency, until incoherent radiation takes over when the wavelength of the 
radiation is much less than electron spacing [Ref. 3: p. 3246]. 

When a single electron passes through a dielectric medium where a 
spatially periodic term is added to the dielectric constant, the result is 
radiation occurring even for electrons which do not exceed the velocity of 
light in the medium and at angles other than the Cerenkov cone angle. The 
non-Cerenkov part of the radiation is attributed to transition radiation. If the 
electron velocity were lower so that v/c were close to but less than unity, the 
peak in I would be pushed to the left toward 0° of 9 such that cos 9 C = v/c 
would be larger than 1. But the tails of the diffraction function I could extend 
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into the physical range 1 = cos 9 — - 1, and this would be called transition 
radiation and be ascribed to the passage of the electrons through the 
boundaries of the gas cell. [Ref. 2: p. 1535] 

F. PROFILE DETERMINATION OF A ELECTRON BEAM BUNCH 

1. Effect of The Form Factor 

The bunch form factors are Fourier Transforms F(k) which differ 
depending on the charge distribution within a bunch. The form factor is 
identically one for a point charge, and for a finite charge distribution F(k) = 
1 for k = 0. If all the charge within a bunch is of the same sign F(k) must 
fall off as a function of k near the origin. Here we consider only line charge 
distribution and consequently the form factors are functions of the product of 
k z = kcos 0 and appropriate lengths which are characteristic of the charge 
distribution of a bunch. 

Thus the behavior of the form factor F over the range of emission angles 
is dependant on the ratio of the charge bunch lengths to the wavelength of 
the emitted radiation. For a small value of these ratios, the value of F is 
essentially unity for all angles. Then the emitted radiation pattern is 
modulated only by the geometrical envelope function as described in Ref. 5. 

However for larger values of the ratios, the value of F over the range of 
emission angles can be significantly less than one. In this case the dependence 
of the form factor on angle is dramatic since the cosine term goes to zero at 
ninety degrees causing a form factor of unity. If, at same time the emission 
length ratio is small, a condition could be achieved where the majority of the 
Cerenkov power is radiated in a direction normal to the electron beam path. 
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A complete angular map of Cerenkov radiation patterns may be used to 
characterize the properties of relativistic beam pulses, and may enable one to 
determine the beam pulse charge distribution from measured radiated patterns 
[Ref. 4: p. 1996]. 

2 Bunch Shapes and Corresponding Form Factors 

Figure 2.3 diagrams the beam pulse shapes and presents the 
corresponding form factor expressions. [Ref. S: pp. 16,17] 



Level function 



(a) 



V :) 




sin k 7 d 

F(k z ) = 

k 2 d 



Figure 2.3 Electron Beam Charge Bunch Shapes and 
Corresponding Form Factors. 
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Triangular function. 



(b) 




Trapezoidal function 




It = top level length of trapezoidal, 
lb = bottom level length of trapezoidal. 

Figure 2.3 (cont.) Electron Beam Charge Bunch Shapes and 

Corresponding Form Factors. 
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Rounded function 



(d) 




F(k 2 ) 



8 

k!(b+d)(d-b) 2 



k z 

2 sin (d+b) - sin k z b - sin k z d 

2 



Gaussian function 




-( kzb ] 2 

F(k z ) = e ' where b = bunch size. 



Figure 2.3 (cont.) Electron Beam Charge Bunch Shapes and 

Corresponding Form Factors. 
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Level Plus Ripple function 




r /, n, sin k z d , R sin(k z + ko)d . sin(k z + ko)d 

t(kz) = + T"l “(k‘ a - k 0 )d + (k z + k 0 )d 

where R = B/A ko = 2x/A 



Multiple Hump function 



(g) 




F(k z ) 



2s in k z d 

zzr 



sin k z d 


2N-3 1 




sin k z d| 

_L 


[ 2N-5] 


2FT J 




k^ 


N=2 kzd 



N=3 



where N = number of humps. 



Figure 2.3 (cont.) Electron Beam Charge Bunch Shapes and 

Corresponding Form Factors. 
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III. CALCULATION AND RESULTS 



A. PARAMETERS FOR CALCULATION 

Professor J. R. Neighbours and others of the Department of Physics, 
Naval Postgraduate School, Monterey, California, developed a computer 
program to calculate and plot angular maps of the radiated energy per unit 
solid angle for a given form factor. The program originally was designed to 
run on a Tektronix Computer and considered the following one dimensional 
(i.e. in Z only) beam charge distributions: Gaussian, level (i.e. boxcar), level 
plus a sinusoidal, ripple combination, and a double hump. 

This work was aimed to modify the program to run on the IBM 370 
mainframe computer and plot the calculated power per solid angle, radiated at 
the frequency v, using the DISSPLA V10.5. The results are limited to the 
three form factors for the rectangular, Gaussian and trapezoidal charge 
distributions. Calculations for other charge distributions are left for future 
work. 

To calculate the total coherent power per unit solid angle, radiated at the 
frequency 1 / by a periodic charged particle beam in traveling a finite distance 
L at constant velocity, a set of parameters were chosen close to those 
available in an experiment. 

(1) Index of refraction, n = 1.000268. 

(2) Speed of light in vacuum, Co = 2.997925 * 10 8 m. 

(3) Fundamental frequency, vo = 2.85 GHz. 

(4) Bunch charge, B = 1.5 * 10' 12 Coulomb. 
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These values give an n P product of n /3 = 0.98, which is 2 % below the 
infinite path length threshold value of unity. By choosing different radiation 
harmonics and varying the path length the radiation pattern above ( 77 = 12 ) 
and below ( 77 = 37 ) the threshold for 2.555 Mev electron bunches, could be 
investigated. The calculations were carried out at constant rj for the various 
charge distribution as the size of the distribution was varied. 

B. GAUSSIAN FUNCTION 

Fig. 3.1 through Fig. 3.8 show the radiation intensity pattern of Gaussian 
line charge distributions in sub-Cerenkov radiation (n/3 = 0.98 ) for two 
different value of 77 . Fig. 3.1 through Fig. 3.4 are the results of using 77 = 12, 
with N = 3, varying the bunch length from 0.24 cm to 2.4 cm. The rest of 
them are for the case of using 77 = 37, and likewise, varying the bunch length. 
BL is the symbol used for bunchlength. 

Figures 3 . 1-3.4 show the effect on the emitted radiation pattern as the 
bunch length is increased while the beam path parameter 77 is held constant at 
77 = 12 . With a bunch length of 0.24 cm the radiation pattern is essentially 
that of a periodic set of point charges. As the bunch length increases, a hump 
first becomes evident in Fig. 3.2 when the wavelength of the emitted radiation 
begins to become comparable to the size of the bunch. Further increase in 
bunch size shows an increasing dominance of the form factor until finally the 
radiation pattern is almost completely dependant on the form factor. In Fig. 
3.4 the envelope of the radiation pattern is the form factor for the Gaussian 
charge distribution except for a slight skewing caused by the geometrical 
envelope function. 
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Fig. 3.1 Gaussian Function with BL/A = 0.069, 77 = 12 
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Fig. 3.2 Gaussian Function with BL/A = 0.343, 7 / = 12 
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Fig. 3.3 Gaussian Function with BL/A = 0.503, 77 = 12 
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Fig. 3.4 Gaussian Function with BL/A = 0.686, 77 = 12 
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Fig. 3.5 Gaussian Function with BL/A = 0.069, 77 = 37 
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Fig. 3.6 Gaussian Function with BL/A = 0.343, rj = 37 
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Fig. 3.7 Gaussian Function with BL/A = 0.503, 77 = 37 
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Fig. 3.8 Gaussian Function with BL/A = 0.686, 77 = 37 
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Similar remarks can be made for the longer finite path results shown in 
Fig. 3.5 - 3.8 for 77 = 37. In this set of figures, because this value of 7? causes 
the radiation to be below threshold, the radiation patterns are wildly 
oscillating - almost chaotic. However exactly the same thing happens to the 
envelope of the radiation. Fig. 3.8, calculated for a bunch length of 2.4 cm, as 
is Fig. 3.4, also shows the radiation envelope to be the form factor for the 
Gaussian charge distribution. 

C. LEVEL FUNCTION 

Fig. 3.9 through Fig. 3.22 show the radiation intensity lobe development 
as the level function charge distribution length varies for the same two values 
of 77 . Fig. 3.9 through Fig. 3.15 are calculated for 77 = 12, the remainder 
through Fig. 3.22 for 77 = 37. LFL is the symbol for the level function length. 

In Fig. 3.9, where the level function length is much smaller than the 
radiation length, the radiation pattern is almost the same shape as Fig. 3.1 
which is for the Gaussian charge distribution. Proceeding through to Fig. 3.13, 
as the length of the charge distribution is increased the form factor becomes 
prominent. Fig. 3.13 is of interest in that it also shows considerable radiation 
at small angles. This result occurs because the form factor for the level 
function has secondary maximum which at low angles combine with the 
relatively large value of the diffraction function to give substantial radiation. 

Further increases in level function length give an increased role to the 
form factor until finally in Fig. 3.15, it has become dominant. Here the 
envelope of the radiation is seen to be the (skewed) form factor of the level 
function. The small lone peak at approximately ten degrees occurs at the 
maximum of the geometrical envelope function. 



34 




Fig. 3.9 Level Function with LFL/A = 0.285, 77 = 12 
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Fig. 3.10 Level Function with LFL/A = 0.855, rj = 12 
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Fig. 3.11 Level Function with LFL/A = 1.28, rj = 12 
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Fig. 3.12 Level Function with LFL/A = 1.71, rj = 12 
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Fig. 3.13 Level Function with LFL/A = 1.71, rj = 12 
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Fig. 3.14 Level Function with LF1/A = 1.71, 77 = 24 
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Fig. 3.15 Level Function with LFL/A = 2.2S, 77 = 12 
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Fig. 3.16 Level Function with LFL/A = 0.29, 77 = 37 
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Fig. 3.17 Level Function with LFL/A = 0.85, 77 = 37 
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Fig. 3.18 Level Function with LFL/A = 1.28, 77 = 37 
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Fig. 3.19 Level Function with LFL/A = 1.71, 7/ = 37 
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Fig. 3.20 Level Function with LFL/A = 1.71, 77 = 37 
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Fig. 3.21 Level Function with LFL/A = 1.71, rj = 74 
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Fig. 3.22 Level Function with LFL/A = 2.28, 77 = 37 
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Fig. 3.16 to Fig. 3.22 drawn for 77 = 37 show the increasing dominance of 
the form factor below threshold. The pattern oscillates very rapidly, but 
otherwise the results are as above. In fact the rapid oscillations make the form 
factor more evident in this series of Figures. 

D. TRAPEZOIDAL FUNCTION 

Both the Gaussian and level function distributions have a single 
characteristic length; the trapezoid is more complicated in having two, the 
length of the base, lb, and the length of the top, It- In the fourier transform 
these appear as ( = ( lb + It )/2 and 1 = ( lb - It )/2. 

Fig. 3.23 through 3.32 show radiation patterns for a trapezoidal charge 
distribution for different values of top and base lengths. All except the last 
three, Fig. 3.30 through Fig. 3.32 are for t] = 12. TFTL is the symbol for the 
top level length of the trapezoidal function and TFBL, for the base level 
length of trapezoidal function. 

As mentioned earlier, the form factor is sensitive to the ratios of the 
characteristic lengths to the wavelength of the emitted radiation. Here the 
appropriate ratios are iff] and 1/A. In Fig. 3.23 both these ratios are less than 
one so that the radiation pattern is characteristic of a point charge 
distribution. Fig. 3.24 drawn for (/ X 1.2 and 1/A = 0.29 shows a different 
pattern with the form factor already becoming quite evident. 

Fig. 3.25 through Fig. 3.27, calculated for a constant slightly larger ratio 
of ifX = 1.71, and increasing values of 1/A are all similar since the radiation 
wavelength is too large to be sensitive to the smaller details of the charge 
distribution. The same comment applies also to Fig. 3.28 for which If X is still 
larger but 1/A is less than one. 
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Fig. 3.23 Trapezoidal Function with = 0.60, 1/A = 0.20 
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Fig. 3.24 Trapezoidal Function with ij A = 1.20, 1/A = 0.29 



51 



W (Walt) 



o 

* 



20.0 



18.0 



16.0 



14.0 - 



12.0 - 



10.0 •- 



8.0 



6.0 



1 

i 

i 

L 



! / 

*T 



1.0 



0.0 



! I 



I 



: i 



Si 



il 



| ;| i) |! 
i , 

* i 









ihln i i 

; 




1 



« ,, I! 

1 ! 



sub-Cerenkov 
E = 2.555 MeV 

T) = 12 

L = 42 cm 
A = 3.50 cm 
TFTL = 5.0 cm 
TFBL = 7.0 cm 



!iS 



! !i 



!i 'I '! 



il !! I 

f h ! i ! 1 1 i 






20 40 60 80 100 120 140 160 180 

Angie (6, degrees) 



Fig. 3.2') Trapezoidal Function with f/X = 1.71, 1/A = 0.29 
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Fig. 3.2C Trapezoidal Function with (/X = 1.71, 1/A = 0.57 
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Fig. 3.27 Trapezoidal Function with (/X = 1.71, 1/A = 0.S6 
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Fig. 3.2S Trapezoidal Function with tj\ = 2.00, 1/A = 0.56 



oo 




Fig. 3.29 Trapezoidal Function with C / A = 1.29, 1/A = 1.27 
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Fig. 3.30 Trapezoidal Function with (j A = 1.20, 1/A = 0.29 




Fig. 3.31 Trapezoidal Function with C/\ = 1.71. 1/A = 0.29 
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Fig. 3.32 Trapezoidal Function with (/\ = 1.71, 1/A = 0.57 
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The last figure of the 77 = 12 series is drawn for both (/X and 1/A greater 
than one. Consequently one would expect the radiation envelope to display 
both terms in the form factors for a trapezoidal charge distribution. That it 
apparently does not may be due to the fact that the lengths chosen in order 
to achieve these make for a charge distribution which is almost triangular. 

Figures 3.30-3.32 are calculated for the below threshold condition of rj — 
37. All have 1 / A greater than one and 1/A less than one. As with other figures 
with this value of 77 , the envelope is more evident. 

When both values of t/X and 1/A increase the peak values decrease mainly 
in the region of angles between 0° and approximately 55° and the radiation 
intensity is enhanced around 90°. It is difficult to determine the development 
of the radiation pattern with respect to if X and 1/A because both values can 
change the radiation pattern at same time. To analyze the radiation pattern in 
another way the angle where zeros happen can be investigated by calculation. 
Let's set: 





m 7T 



(3.1) 



where k z = k cos 0, and m,n are integers. Substituting for gives 



cos $ n 




cos $ m = ± m -A 



A 



(3.2) 



If as in most cases considered 




1, these are no solutions for 6 m . 
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Thus only the solutions for 9 a are considered below. A example is Fig. 3.25 



c 

where — ^ — — 1.71, i — 6 cm, A — — 3.5 cm. 

cos 0 = ± n = ± n = ± n * 0.58 (3.3) 

The solutions are for 9 n are 



0i , 54.5° and 9. i = 125.5°. 



(3.4) 



which agree with the Figure. These values also apply to Fig. 3.26, Fig. 3.27, 

f 

Fig. 3.31, and Fig. 3.32 which have the same ratio of — r-. 

In Fig. 3.29 where -j- 1.29, 1/A = 1.27, 



cos 9 n = ± n 



1.29 



COS 9m 



± m 



1.27 



(3.5) 



The solutions are = 39.2°, = 129.2° for 0 n and 6\ = 38°, 9. i = 128°. In 

This case the zeros for 9 m and 9 n are so close together as to be 
indistinguishable in the figure. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The object of this work is to assess the effects of the charge distribution 
of a bunch on the sub-Cerenkov radiation from a charged particle beam, 
specifically investigating the short beam path region at energies above and 
below the threshold. It contains many calculated radiation patterns for the 
three kinds of charge distributions. 

With the charge distribution of Gaussian function, the radiation intensity 
envelope near 90° represents the square of the bunch form factor F(k) when 

the bunch length is larger than a half of the wavelength. When the bunch 

length is about same as the wavelength or larger than that of the level charge 
distribution the radiation pattern around 90° is also dominated by the square 
of the bunch form factor F(k). 

It is difficult to determine the characteristics of the radiation pattern on 
the trapezoidal function by the ratio of bunch length to wavelength because 

the level top length and level base length can change the shape of radiation 

pattern at same time. However the radiation pattern of the trapezoidal 
function can be partially characterized by analyzing the zeros of the curves. 

A summary of results is : 

1. When the bunch length is much smaller than the wavelength of the 
emitted radiation, the radiation pattern is characteristic of point 
charge. 

2. When rj is increased the total radiation pattern does not change 
except for minor differences in the frequency of oscillation. 
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3 . 



When the wavelength of the emitted radiation is comparable to or 
larger than the bunch length, a significant amount of radiation 
appears at near 90° to the beam. 

4. The envelope of the 90° radiation is the square of the bunch form 
factor. This envelope is skewed slightly as a result of monotonic character of 
the geometrical envelope function. 

To enhance this preliminary work the following is recommanded: 

1. Investigate the behavior of the square of form factor F(k) only when 
the bunch length varies. This will give the insight into the behavior 
of radiation pattern. 

2. Determine why varying the top length of the trapezoidal pulse shape 
causes significant enhancement of the forward lobe envelope. 

3. Determine the effects of changes in beam energy, pulse width, pulse 
frequency, and path length. 

4. Confirm the results found by mapping the microwave radiation 
measured at a comparable particle accelerator. 

The examples presented here are idealized and there has been no attempt to 
take into consideration such factors as beam pulse erosion as it transverses the 
atmosphere and the ground plane reflection in interference[Ref. 4: p. 1996] . But 
the complete angular map of the envelope of the radiation pattern could 
provide a measure of F(k z ) and an excellent reference to establish a method of 
determining the beam pulse charge distribution based on sub-Cerenkov 
radiation patterns. 
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APPENDIX A 



PROGRAM FOR CALCULATION OF RADIATED POWER 



********************#**##***#***•##••#-*•#.*•*-*#■*•#-*•#•**•*#■.*•*-*•**###.*##.*#*.#■#*#.#■#* 

* THIS PROGRAM IS TO CALCULATE THE POWER PER UNIT SOLID ANGLE RADIA- * 

* TED AT THE FREQUENCY NU. IT WAS ORIGINALLY DEVELOPED BY PROFESSOR * 

* JOHN R. NEIGHBOURS AND OTHERS FOR TECTRONIX COMPUTER. THIS AUTHOR * 

* MODIFIED TO RUN ON IBM 370 MAINFRAME COMPUTER IN DEC. 1987. * 

* THE THEORY AND ALL EQUATION ARE WELL GIVEN IN REF. 5. * 

*********************■**■*■#•*•#*###*#.**##*#.##•*.#■**.*#.###••*#•#*###.*#•#•#•#*•#.*.##.#.#.# 



C 

C 

c 



c 

c 

c 



PROGRAM JNN 

COMMON / VA/HH , HU , SDXX /DOUP A/L , ETA , TMAX , DT , E , N , B , LL , RR , NN , RFTL , 
(- RFBL , TFTL , TFBL , TRIBL , NBETA 

C , C0 , MU0 , NU0 , Q , P I , E0 , E , NBETA , 

TMAX , DT , NU , L , ETA , N , B , LL , RR , NN 
UMAX , KK , CONST , WMAX , THETA , W , T , KKZ , H , DENOM1 , 
DENOM2 , K0 

GAMMA , V , U , I I , R , HA , HB , HC , HD , DD , DEXIN 
RFTL , RFBL , TFTL , TFBL , TRIBL 

DENOM3 , DEN0M4 , DENOM5 , DENOM6 , DEN0M7 , DEN0M8 , 

HE , HF , HG , HI 

H J , HK , HL , Z , L I , L IN , HX , HW , KE , SDXX 
XWMAX , XTHETA , XW , YMAX , SDYY , SDY , JX , J YY , YY , YNUM , X , Y , SDX , 

+ XX , XNUM , RTMAX , RSDXX , RWMAX , RTHETA , RW 

REAL REL , REE , REN , REB , RELL , RERR , RENN , DELTA 
REAL RERFTL , RERFBL , RETFTL , RETFBL , RETRIB , REHU 
INTEGER I , K , KMAX , IMAX , HH , J , JXMAX , JYMAX , IP A , P , HU 
CHARACTER* 1 P ARACH , EXPPTS , EXPTCH , SAVEQ , PPRE , CHANGE , ANS , AXCH , 

+ ANS2 , ANS3 

CHARACTER* 1 3 NAMEPF , NAMEWF , NAMEEX 
DIMENSION XTHETA( 1:100) 

DIMENSION XW( 1 : 100) 

DIMENSION THETA( 0 : 4500 ) 

DIMENSION W( 0 : 4500 ) 

DIMENSION H( 1 : 7 ) 



DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
REAL 



* ASSIGNMENT OF CONSTANTS 



C=2 . 9971 22D08 
C0=2 . 997925D08 
MU0= 1 . 256643706 1 D-06 

********** rpo VARY CHARGE ************ 



PI=3. 141592653 5897932D0 
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E0 = 0 . 5 1 1 7D0 



5 PRINT* , ’ 

PRINT*, ’ ENTER BUNCH CHARGE .( DEFAULT= 1.5E-12) ’ 

READ*, Q 

IF( Q . GT . 0 ) GO TO 45 
Q= 1 . 5D- 1 2 

45 PRINT*, ’ ENTER FUNDAMENTAL FREQUNCY NU0 . ( DEFAULT= 2.85E9)’ 

READ*, NU0 

IF( NU0 . GT . 0 ) GO TO 55 
NU0=2 . 85D09 
55 CONTINUE 

C 



C 

C 

C 



PRINT* , * THE INDEX OF REFRACTION IN THE AIR IS 1 . 000268( DEFAULT ). 1 0 
PRINT* , ’ DO YOU WISH TO CHANGE? (Y/N) ’ 

READ 1 ( 1 A ) * , ANS2 
IF( ANS2 . EQ . * Y * ) THEN 

PRINT* , ’ENTER THEINDEX OF REFRACTION. 

READ*, DEXIN 
ELSE 

DEXIN= 1 . 000268D0 
END IF 
C=C0/DEXIN 



* * * * * * * 



DECISION OF PARAMETER FILE******* 



PRINT * 

PRINT *, ’ 

PRINT *, ’ 

PRINT *, ’ 

PRINT * 

PRINT *, ’ 

READ ’ ( 1 A) ’ , 
IF( PARACH . EQ , 



1: CREATE NEW PARAMETER FILE 1 
2: USE OLD FILE’ 

5: EXIT’ 

ENTER YOUR CHOICE : * 

PARACH 

* 5 * ) GOTO 700 



IF( PARACH . EQ . ’2’ ) CALL INAMPF 
IF( PARACH . EQ . • 1 * ) CALL MAKEPF 



C 

***************************************************** 



* READ IN ALL PARAMETERS FROM EXTERNAL * 

* FILE/PRINTOUT PARAMETERS 
***************************************************** 



C 



99 


WRITE( 6 , 1 00) 




1 00 


FORMAT( III ' THE PARAMETERS YOU ENTERED 
REWIND 1 


ARE: ’ ) 


1 20 


READ( 1 , * ) HH 




151 


WRITE( 6,135) HH 




155 


FORMAT( / ’FOURIER TRANSFORM CHOICE : 
IF (HH.NE. 1 ) GOTO 1 37 
WRITE (6,136) 


’ , 18) 


156 


FORMAT(/ ,37X, ’GAUSSIAN FUNCTION ’ ) 




1 57 


IF (HH.NE. 2) GOTO 139 
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1 38 

139 

140 

141 

142 

143 

144 

1 45 

146 

147 

148 

C 

149 

9151 

1 51 

154 

155 

1 64 

165 

174 

175 

C 

Cl 84 

184 

Cl 85 

185 

1 94 

195 

C 

C 

C 

204 

205 

C 

206 

C 

214 



WRITE (6,138) 

FORMAT (/ ,37X, ’LEVEL FUNCTION ’ ) 

IF (HH.NE.3) GOTO 141 
WRITE (6,140) 

FORMAT (/’SUM OF LEVEL AND RIPPLE FUNCTIONS’ ) 

IF (HH.NE.4) GOTO 143 
WRITE (6,142) 

FORMAT( / ’MULTIPLE HUMP FUNCTION’ ) 

IF (HH.NE.5) GOTO 145 
WRITE (6,144) 

FORMAT (/’ROUNDED FUNCTION' ) 

IF (HH.NE.6) GOTO 147 
WRITE (6,146) 

FORMAT (/’TRAPEZOIDAL FUNCTION’ ) 

IF (HH.NE.7) GOTO 149 
WRITE (6,148) 

FORMAT (/’TRIANGLE FUNCTION’ ) 

READ( 1 , * ) ETA 
WRITE(6, 151 ) ETA 

FORMAT( / ’ ETA ( L/RAMMDA ) : ’.D16.4) 

READ( 1 , * ) TMAX 
WRITE( 6,155) TMAX 

FORMAT ( / ’ MAXIMUM THETA VALUE : ’.D16.4) 

READ( 1 ,*) DT 
WRITE( 6,165) DT 

FORMAT ( / ’ THETA CALCULATION INCREMENT :’,D16.4) 
READ( 1 , * ) SDXX 
WRITE( 6,175) SDXX 

FORMAT ( / ' THETA MARKING INCREMENT : ’.D16.4) 

READ( 1 , * ) E 
READ( 1 ,*) NBETA 
WRITE( 6,185) E 
WRITE( 6,185) NBETA 

FORMAT ( / ’ BEAM ENERGY : \D16.4) 

FORMAT ( / ’ VALUE OF NBETA : ’ , D 1 6 . 4 ) 

READ( 1 , * ) N 
WRITE( 6,195) N 

FORMAT ( / ’ HARMONIC NUMBER : ’.D16.4) 



IF (HH.NE. 1 ) GOTO 206 



READ( 1 ,*) B 

WRITE( 6,205) B 

FORMAT ( / ’ BUNCH LENGTH (M) 

IF (HH.LT.2 ) GOTO 216 
IF (HH.GT.4) GOTO 216 

READ( 1 , * ) LL 
WRITE( 6,215) LL 



’ , D1 6 . 4 ) 
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215 FORMAT ( / ’ LEVEL FUNCTION LENGTH (M) : ' ,D16.4 

C 

216 IF (HH.NE.5) GOTO 9255 
C 

READ( 1 ,*) RR 

224 WRITE( 6 , 225 ) RR 

225 FORMAT ( / ’ FOURIER RATIO CONSTANT : ' ,D16.4) 

C 

C 

READ( 1 ,*) NN 

254 WRITE( 6,255) NN 

255 FORMAT ( / ' NUMBER OF RIPPLES : ' ,D16.4) 

C 

9255 IF (HH.NE.4) GOTO 256 
C 

READ (1 ,*)HU 

9257 WRITE (6,9259) HU 

9259 FORMAT ( / ’ NUMBER OF HUMPS : ’,14) 

C 

256 IF (HH.NE.5) GOTO 241 
C 

READ ( 1 , * ) RFTL 

9258 WRITE (6,258) RFTL 

258 FORMAT (/’ ROUNDED FUNCTION LEVEL TOP LENGTH : 

C 

READ ( 1 , * ) RFBL 
9240 WRITE (6,240) RFBL 

240 FORMAT (/’ ROUNDED FUNCTION BASE LENGTH : 

C 

241 IF (HH.NE.6) GOTO 246 
C 

READ (1 ,*) TFTL 
9245 WRITE (6,245)TFTL 

245 FORMAT (/’ TRAPEZOIDAL FUNCTION TOP LEVEL LENGTH : 
C 

READ ( 1 , * ) TFBL 
9245 WRITE (6,245)TFBL 

245 FORMAT (/’ TRAPEZOIDAL FUNCTION BASE LENGTH : 

C 

246 IF (HH.NE.7) GOTO 249 
C 

READ (1 ,*) TRIBL 
9248 WRITE (6,248)TRIBL 

248 FORMAT ( / ’ TRIANGLE FUNCTION BASE LENGTH : 

249 WRITE(6,250) 

250 FORMAT ( / / ’ DO YOU WISH TO CHANGE ANY OF THESE ?(Y OR 
READ * ( A 1 ) ’ , CHANGE 

C 

IF( CHANGE. EQ. *Y’ ) GO TO 700 
C 

-it****************************** 

* *** BEGIN CALCULATION *** * 



’ ,D1 6.4) 
’ ,D1 6.4) 

’ , D1 6 . 4 ) 
■ ,D1 6.4) 

’ , D1 6 . 4 ) 
N)’ ) 
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#*#-***-»**#*##*•***#####*-*#**#**# 



c 

325 IIMAX=TMAX/DT 
IMAX= INT ( UMAX) 

NU=N*NU0 

IF ( ( HH . GE • 2 ) . AND . ( HH . LE . 4 ) ) THEN 
DD=LL/ 2 . 0D0 
END IF 
C 

KK=2 . 0D0*P I*NU / C 

CONST =MU0*C* ( NU0**2 .0)*(Q**2.0)/ (8.0D0*(PI**2.0) ) 

WMAX=0.0D0 

DO 350 1=0, IMAX 

THETA( I ) =DBLE( I ) *DT 
T=PI*THETA( I ) / 1 80 . 0D0 
KKZ=KK*DCOS( T ) 

C 

IF ( HH.NE. 1 ) GOTO 326 
C 

H( 1 ) =DEXP ( - 1 .0D0*((KKZ*B/2,0D0)**2.0)) 

C 

326 CONTINUE 

IF (HH.LT.2) GOTO 328 
IF (HH.GT.3) GOTO 328 
C 

DENOM1 =KKZ*DD 

IF ( DENOM1 . EQ . 0 ) GOTO 327 

H( 2 ) = ( DSIN ( KKZ*DD ) ) / ( KKZ*DD ) 

GOTO 328 

327 H( 2 ) =0 . 0D0 
C 

328 IF( HH . NE . 3 ) GOTO 333 

C 

KO=PI*NN/DD 
DENOM2= ( KKZ-KO ) *DD 
C 

IF( DEN0M2 . EQ . 0 ) GOTO 329 
C 

HA= DSIN( DENOM2 ) /DEN0M2 
C 

GOTO 330 
C 

329 HA= 0 . 0D0 

330 DENOM3= ( KKZ+KO ) *DD 
C 

IF (DENOM3.EQ.0) GOTO 331 
C 

HB=DSIN( DENOM3 )/DENOM3 
C 

GOTO 352 
C 

331 HB=0 . 0D0 
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332 
C 

333 
C 

C 

C 



9333 

C 

C 

334 
C 

335 
C 



C 

339 

C 



C 

345 

C 

C 



C 

C 

346 



H(3)=H(2)+((RR/2.0D0)*( HA+HB ) ) 

IF( HH . NE . 4 ) GOTO 335 

Z=(KKZ*DD) 

IF( Z . EQ . 0 ) GOTO 334 

H( 4 ) = 1 . 0D0*( ( DSIN( Z ) )/Z ) 

DEN0M4= 1 .0 
DO 9333 P=1 , HU 

KE= ( 2 . 0D0*DBLE( P ) )- 1 . 0D0 

LI =( ( 2 . 0D0*DBLE( HU ) )-KE ) / ( ( 2 . 0D0*DBLE( HU ) )- 1 . 0D0 ) 

LIN=( ( -1 .0D0)**DBLE(P))*(-1 . 0D0 ) 

HX=(LIN*LI*DSIN(Z*LI ) )/(Z*LI) 

HW=LIN*LI 
H( 4 ) =H( 4 )+HX 
DENOM4= DEN0M4+HW 
CONTINUE 

H( 4 ) =H( 4 ) /DEN0M4 
GOTO 335 
H( 4 ) =0 . 0D0 

IF (HH.NE.5) GO TO 339 

HL = ( 8D0 * 1 6D0 ) / ( ( ( RFBL-RFTL ) * *2 . 0 ) * ( KKZ** 3 ) * ( RFBL+RFTL ) ) 
HG=HL*DSIN( ( KKZ / 4 . 0D0 ) * ( RFBL+RFTL ) ) 

HJ = (RFTL*DSIN(KKZ*RFTL/2.0D0) ) / ( KKZ*RFTL/2 . 0D0 ) 

HK= (RFBL*DSIN(KKZ*RFBL/2 . 0D0 ) ) / (KKZ*RFBL/2 . 0D0 ) 

HI = ( 32D0 / ( ( ( ( RFBL-RFTL ) * ( KKZ ) ) **2 . 0 ) * ( RFTL+RFBL ) ) ) * ( H J+HK ) 
H( 5 )=HG-HI 

IF (HH.NE.6) GOTO 345 

DENOM6=KKZ*( TFBL+TFTL ) /4 . 0D0 
DEN0M7=KKZ* ( TFBL-TFTL ) / 4 . 0D0 
HE=DSIN( DENOM6 ) /DENOM6 
HF=DSIN( DEN0M7 ) /DEN0M7 
H( 6 )=( TFBL+TFTL ) /2 . 0D0*HE*HF 

IF( HH . NE . 7 ) GOTO 349 

DENOM5=(KKZ*TRIBL)/4.0D0 

IF (DENOM5.EQ.0) GOTO 346 

H( 7 ) = ( DSIN( DENOM5 ) / ( DENOM5 ) ) **2 . 0 

GOTO 349 

H( 7 ) =0 . 0D0 
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CONTINUE 



C 

349 
C 

C GAMMA=E/E0 

L=ETA*C/NU 

C GAMMA=DSQRT( 1 .0D0/(1 . 0D0-( NBETA/DEXIN ) * *2 ) ) 

C V=C0*( DSQRT( 1 . 0D0-( 1 . 0D0 / ( GAMMA**2 . 0 ) ) ) ) 

V=C0*( NBETA/DEXIN ) 

U=(KK*L/2.0D0)*( (C/V)-DCOS(T) ) 

II=DSIN(U)/U 

R=KK*L*H(HH)*II*DSIN(T) 

W(I)=CONST*(R**2.0) 

C 

C *** CALCULATION OF --MAXIMUM W *** 

C 

272 WMAX=DMAX1 (WMAX , W( I ) ) 

C 

350 CONTINUE 
C 

C ****** DECISION TO SAVE VALUES CALCULATED ****** 

C 

WRITE( 6 , 360 ) 

360 FORMAT ( / ’ DO WANT TO SAVE THE RESULTS IN A FILE? (Y/N) 1 ) 

READ' ( A1 ) ’ , SAVEQ 
C 

IF( SAVEQ. EQ. ’N’ ) GOTO 440 
C 

PRINT *, 'UNDER WHAT FILE NAME? ' 

READ ' ( A1 0 ) ’ , NAMEWF 

OPEN( UNIT = 3 ,FILE=NAMEWF, STATUS= ' NEW’ , FORM= ' FORMATTED ' ) 

DO 370 1=0 , IMAX 

RTHETA=REAL( THETA( I ) ) 

RW=REAL(W( I ) ) 

WRITE( 3 , 380 ) RTHETA , RW 



380 


FORMAT (4X,F8. 3, 


2X.F20.18) 


370 


CONTINUE 
WRITE( 3 , 385 ) 




385 


FORMAT( 7X , '0.000' 
WRITE( 3 , 390 ) WMAX 


) 


390 


FORMAT (2X, 'WMAX = 
WRITE( 3 , 392 ) ETA 


' ,D1 6.4) 


392 


FORMAT (2X, 'ETA = 


’ ,D1 6.4) 




WRITE( 3,394) NBETA 


394 


F0RMAT(2X, 'NBETA 


= ’ , D1 6 . 4 ) 


C 


WRITE( 3 , 396 ) B 




C396 

r* 


F0RMAT(2X, 'BUNCH 


LENGTH = ’ ,D1 6.4) 


C ****** DECISION TO 

r* 


PRINTOUT, RE-RUN, OR EXIT 


440 


WRITE( 6 , 445 ) 




445 


FORMAT(////’ 1 


: PRINT OUT VALUES’,/ 




+ ' 2 


: RUN PROGRAM AGAIN’ ,/ 
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+ ’5: EXIT’ ,/// 

+ ' ENTER CHOICE. ’) 

READ’ ( A1 ) ’ , PPRE 
C 

IF( PPRE . EQ . ’ 1 ’ ) GOTO 660 
IF( PPRE . EQ . ’2 ’ ) GOTO 05 
IF( PPRE .EQ . ’ 5 ’ ) GOTO 700 
GOTO 440 
C 

C *** PRINT OUT VALUES *** 

C 

660 WRITE( 6,670) 

670 FORMAT ( ’ THETA ¥ , /23(’-’)/) 

DO 690 1=0, IMAX 

WRITE( 6 , 680 ) THETA( I ) , W( I ) 

680 FORMAT ( F8 .2 , FI 5 • 4 ) 

690 CONTINUE 
C 

C ********* PROGRAM ENDS ********************************** 

700 END 

C 

C 

c 

*'l***#***'l*****#*»**»****##»**# 

* *** SUBROUTINE INAMPF **» * 

»•»**•**»•***•****•*■»**■»*********'* 

C 

SUBROUTINE INAMPF 
C 

C *** INPUT OLD PARAMETER FILE NAME 
CHARACTER* 10 NAMEPF 
WRITE( 6,900) 

900 FORMAT( / ’ENTER THE NAME OF THE PARAMETER FILE YOU WISH TO’, 
+ ’ USE( MAX 10 CHAR) ’ ) 

READ ’ ( A 1 0 ) ’ , NAMEPF 
WRITE( 6,91 1 ) 

911 FORMAT( / ’TEST LINE ONE’ ) 

OPEN( UNIT= 1 , FILE=NAMEPF , STATUS= ’ OLD ’ , FORM= ' FORMATTED ’ , 

+ BLANK= ’ NULL ’ ) 

WRITE( 6,912) 

912 FORMAT( / ’TEST LINE TWO’ ) 

RETURN 

END 

C 

***************************** 

* *** SUBROUTINE INAMEX *** * 

***************************** 

C 

SUBROUTINE INAMEX 
C 

C *** INPUT OLD EXPERIMENTAL POINTS FILE NAME *** 

CHARACTER* 1 0 NAMEEX 
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o o 



1 000 



C 

C 

# # * * # 
* * * # 
* # * * * 
c 



c 



1 300 



1 320 



1330 



1340 

C 

1356 
C 

1357 



C 



WRITE( 6 , 1 000) 

FORMAT( / ’ENTER NAME OF EXPERIMENTAL POINTS FILE YOU WISH’, 
’TO USE( MAX 10 CHAR) ’ ) 

READ * ( A1 0 ) * , NAMEEX 

OPEN ( UNIT = 2 , FILE = NAMEEX , STATUS = ’ OLD * , FORM= ’ FORMATTED ’ , 
BLANK= ’NULL * ) 

RETURN 

END 



*##*###*****#*■#*#***•#*■#* 

SUBROUTINE MAKEPF *** * 
**»-*********#*********-** 



SUBROUTINE MAKEPF 



*** INPUT OF PARAMETERS INTO NEW FILE (UNIT=1)*** 



COMMON / VA/HH , HU , SDXX /DOUPA/L , ETA , TMAX , DT , E , N , B , LL , RR , NN , RFTL 
+ RFBL , TFTL , TFBL , TRIBL , NBETA 

DOUBLE PRECISION L , ETA , TMAX , DT , E , N , B , LL , SDXX , NBETA , 

+ RR , NN , RFTL , RFBL , TFTL , TFBL , TRIBL 

INTEGER HH.HU 
CHARACTER* 10 NAMEPF 
WRITE( 6 , 1 300 ) 

FORMAT( / ’ENTER THE NAME FOR THE NEW PARAMETER FILE.’, 

+ ’(MAX 10 CHAR) ’ ) 

READ ’ ( A1 0 ) ’ , NAMEPF 

OPEN ( UNIT= 1 , FILE=NAMEPF , STATUS= ’ NEW’ , FORM= ’ FORMATTED ’ , 

+ ACCESS= ’ SEQUENTIAL ’ , BLANK= ’ NULL ’ ) 

PRINT*, ’ 

WRITE( 6 , 1 320) 

FORMAT( /////’ NOTE: ALL THE PARAMETERS MUST BE ENTERED’, 

+ ’ ALTHOUGH SOME MAY NOT BE USED’) 

WRITE( 6,1330) 



FORMATC / ’ 


1 : 


GAUSSIAN FUNCTION’/’ 


2: LEVEL FUNCTION 


+ 


r 


3: 


SUM OF LEVEL AND RIPPLE 


FUNCTIONS’ , 


+ 


r 


4 : 


MULTIPLE HUMP FUNCTION’ 


* 


+ 


r 


5: 


ROUNDED FUNCTION’ , 




+ 


r 


6: 


TRAPEZOID FUNCTION’, 




+ 


r 


7: 


TRIANGLE FUNCTION’ ) 





WRITE( 6 , 1 340 ) 

FORMAT( / ’ENTER THE FOURIERTRANSFORM CHOICE:’) 
READ*, HH 



IF( HH . LE . 7 ) GO TO 1359 
WRITE( 6 , 1 357 ) 

FORMATC FOURIER TRANSFORM CHOICE NOT DEFINED.’///, 

+ ’ ENTER NEW CHOICE: ’ ) 

READ*, HH 
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GO TO 1356 



C 

1359 



WRITE( 1 , *) HH 
WRITE( 6 , 1 370) 

1370 FORMAT( / ’ENTER ETA (L/RAMMDA). ( DEFAULT = 12.0D-2 ) ’ ) 

READ*, ETA 

IF( ETA . GT . 0 . 0 ) GOTO 1385 
ETA= 1 2 . 0D0 

1 385 WRITE( 1 ,*) ETA 
WRITE( 6,1400) 

1400 FORMAT( / ’ENTER THE MAXIMUM VALUE OF THETA TO BE’, 

+ ’ CALCULATED(DEG) . ’ / ’ ( DEFAULT = 45.0 ’) 

READ*, TMAX 

IF ( TMAX . GT . 0 . 0 ) GOTO 1415 
TMAX=45 • 0D0 
1415 WRITE(1 ,*) TMAX 
WRITE( 6,1430) 

1430 FORMAT(/ ’ENTER THE ANGLE CALCULATION INCREMENT( DEG) ’ , 

+ ’ ( DEFAULT = 0 . 1 ) ’ ) 

READ*, DT 

IF( DT . GT . 0 . 0 ) GOTO 1445 
DT=0 . 1 0D0 

1445 WRITE( 1 ,*) DT 
WRITE( 6 , 1 460 ) 

1460 FORMAT( / ’ENTER THE ANGLE MARKING INCREMENT (DEG) ’ , 

+ ’ ( DEFAULT = 1 5 • 0 ) ’ ) 

READ*, SDXX 

IF( SDXX . GT . 0 . 0 ) GOTO 1475 
SDXX=1 5.0D0 
1 475 WRITE( 1 ,*) SDXX 
WRITE( 6,1490) 

Cl 490 FORMAT( / ’ENTER THE BEAM ENERGY (MEV) ( DEFAULT= 100.0) ’ 
1490 FORMAT( / ’ENTER THE VALUE OF NBETA ( DEFAULT = 0 . 98 ) ’ ) 

C READ* , E 

READ*, NBETA 

C IF( E . GT . 0 . 0 ) GOTO 1505 

IF( NBETA. GT. 0.0) GOTO 1505 
C E= 1 00 . 0D0 



E=0 . 98D0 

Cl 505 VRITE( 1 ,*) E 
1 505 WRITE( 1 , * ) NBETA 
VRITE( 6 , 1 520 ) 

1520 FORMAT( / ’ENTER THE HARMONIC NUMBER ( DEFAULT= 1.0) ’ ) 

READ*, N 

IF( N . GT . 0 . 0 ) GOTO 1535 
N= 1 . 0D0 

1 535 WRITE( 1 ,*) N 
C 

IF (HH.NE. 1 ) GOTO 1 575 
C 

WRITE( 6,1550) 

1550 FORMAT( / ’ENTER THE BUNCH LENGTH (M).( DEFAULT = 0.24D-2 ) 



) 
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c 



READ*, B 



C 

1565 

C 

1575 

C 

1 580 



C 

C 

1595 

C 

1 605 
C 

1610 

C 

C 

1 625 
C 

1655 

C 

1 640 

C 

C 

1655 

C 

1665 

C 

1 670 

C 

C 

1 680 
C 



IF( B . GT . 0 . 0 ) GOTO 1565 



B=0 . 24D0 
WRITE( 1 ,*) B 

IF ( (HH.LT.2) .OR. (HH.GT.4) ) GOTO 1605 
WRITE( 6 , 1 580 ) 

FORMAT( / ’ENTER THE LEVEL FUNCTION LENGTH ( M) . ’ , 

+ ’ ( DEFAULT=0 .01) ’ ) 

READ*, LL 

IF( LL . GT . 0 . 0 ) GOTO 1595 

LL=0 ♦ 0 1 D0 
TJ RITE( 1 ,*) LL 

IF (HH.NE.5) GOTO 1655 

WRITE( 6,1610) 

FORMAT( //// ’ENTER THE FOURIER RATIO CONSTANT ( DEFAULT= 1 . 0 ) . ’ ) 
READ*, RR 

IF(RR.GT.0.0) GOTO 1625 



RR= 1 . 0D0 
WRITE( 1 ,*) RR 

IF (HH.NE.5) GOTO 1665 

WRITE( 6 , 1 640 ) 

FORMAT(/ ’ENTER THE NUMBER OF RIPPLES ( DEFAULT=1 . 0 ) ’ ) 

READ*, NN 

IF( NN . GT . 0 . 0 ) GOTO 1655 

NN = 1 .0000 
WRITE( 1 ,*) NN 

IF (HH.NE.4) GOTO 1700 

WRITE( 6 , 1 670 ) 

FORMAT( / ’ENTER NUMBER OF HUMPS ( DEFAULT=2 ) ’ ) 

READ*, HU 

IF( HU . GT . 2 ) GOTO 1680 



HU = 2 

VRITE( 1 ,*)HU 
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1700 IF (HH.NE.5) GOTO 1800 
C 

WRITE (6,1710) 

1710 FORMAT (/ ’ENTER THE ROUNDED FUNCTION LEVEL TOP LENGTH (CM)’ 
+ / ’ (DEFAULT = 1.0) ’ ) 

READ*, RFTL 
C 

IF ( RFTL. GT. 0.0) GOTO 1730 
C 

RFTL = 1 . 0D0 
1730 WRITE ( 1 ,*) RFTL 
WRITE (6,1760) 

1760 FORMAT (/’ENTER THE ROUNDED FUNCTION BASE LENGTH (CM).’/, 

+ ’ (DEFAULT =2.0) ’ ) 

READ*, RFBL 
C 

IF ( RFBL. GT. 0.0) GOTO 1780 
C 

RFBL = 1 . 0D0 
1780 WRITE (1,*) RFBL 
C 

1800 IF (HH.NE.6) GOTO 1900 
C 

WRITE (6,1810) 

1810 FORMAT (/'ENTER THE TRAPEZOIDAL FUNCTION LEVEL TOP’, 

+ ’LENGTH (CM).’/’ (DEFAULT = 1.0)’ ) 

READ*, TFTL 
C 

IF ( TFTL. GT. 0.0) GOTO 1830 
C 

TFTL = 1 .0D0 
1830 WRITE ( 1 ,*) TFTL 
WRITE (6,1860) 

1860 FORMAT (/’ENTER THE TRAPEZOIDAL FUNCTION BASE LENGTH (CM)’ 
+ /’ (DEFAULT =2.0) ’ ) 

READ*, TFBL 
C 

IF ( TFBL. GT. 0.0) GOTO 1880 
C 

TFBL = 1.0 

1880 WRITE (1,*) TFBL 
C 

1900 IF (HH.NE.7) GOTO 1990 
C 

WRITE (6,1910) 

1910 FORMAT (/’ENTER TRAIANGLE FUNCTION BASE LENGTH (CM).’//, 

+ ’ (DEFAULT = 1.0) ’ ) 

READ*, TRIBL 
C 

IF ( TRIBL. GT. 0.0) GOTO 1930 
C 

TRIBL = 1 . 0D0 
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1930 


WRITE ( 1 ,*) TRIBL 
WRITE (6,1941) 




1941 


FORMAT(/ ’TEST LINE 


THREE 


C 






1 990 


ENDFILE 1 




C 


WRITE( 6 , 1 991 ) 




1991 


FORMAT (/ 'TEST LINE 

RETURN 

END 


FOUR’ 


C 


-■* 
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APPENDIX B 



PROGRAM FOR PLOTTING THE RADIATION PATTERN 



****##*#****#***********#************************************-***-**-**** 

* THIS PROGRAM IS TO PLOT THE POWER PER UNIT SOLID ANGLE RADIATED * 

* AT THE FREQUENCY NU. THE THEORY AND ALL EQUATIONS ARE DESCRIBED 

* IN REF. 3. IT ASSUMES THETA IN COL. 1 (DEGREES) AND POWER IN COL. 2 * 

* (WATT) RESPECTIVELY. IT SUPPORTS A SINGLE PLOT OF MULTIPLE CURVES * 

* AND USES THE COMPRS/DISSPOP OPTION OF DISSPLA. * 

***##*****************#****#**#*#********#****■*■*■#**•#**■*■***■**•*#*■#•*•#*##••# 

C 

REAL YX , XMIN , XMAX , YMIN . YMAX , X , Y , DX , AX , AY , XPHY , YPHY 
INTEGER IOUT , ITICK, ITICKX , ITICKY, IYVAR , ICT , INN, NB , NETA , NF , 

+ NBL , LF , NMOD , NRAM , TLL , BL , NL 

CHARACTERS YES.ANS 
CHARACTER*4 BLNK( 2 ) 

CHARACTER*60 IMESS , IPHI , IRAM , IMOD , INB , IBL , IETA , LFL , INF , ITLL , 

+ INL 

C 

DIMENSION X( 2000 ) , Y( 2000 ) 

DIMENSION YX( 1 : 5000 ,1:2) 

C 

DATA YES/ * Y ’ / 

DATA BLNK/ ’ ' , ' ' / 

EPS= . 00 1 
C 

C INITIALISE PLOTTER 

C 

100 CONTINUE 

9400 FORMAT(AI) 

9401 FORMAT ( I 1 ) 

9402 FORMAT ( 12 ) 

WRITE( 4 , 9403 ) 

9403 FORMAT ( ’ DO YOU WANT:’,/, 

X ’ 1 . COMPRES’ ,/ , 

X ’ 2. TEK6 1 8 ’ ,/ , 

X ’ ENTER NUMBER (II)’,/) 

READ( 4 , * ) IOUT 
IF( IOUT.EQ. 1 ) THEN 
CALL COMPRS 

ELSEIF( IOUT.EQ. 2) THEN 
CALL TEK618 
ENDIF 

AX = 4.6810 
AY = 6.7416 
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CALL HWROT( ’MOVIE' ) 

CALL PAGE( 8.5,1 1.0) 

XPHY=2 . 0840 
YPHY=2 .7953 
C 

C DEFINE: LOVER/UPPER CASE STD AND GREEK SHIFT LABELS (+,!,?,*). 
C SHIFT FOR INSTRUCTION SET(&), AND THE TYPE FACE. 

C 

CALL MX1 ALF( ’STANDARD’ , ’ ! ’ ) 

CALL MX2ALF( ’L/CSTD’ , ' + ’ ) 

CALL MX3ALF( ’GREEK’ , ’?’ ) 

CALL MX4ALF( ’L/CGREEK ’ , ) 

CALL MX5 ALF( ’ INSTRU ’,’&’) 

CALL SWISSL 

CALL SHDCHR( 90 ., 1 , 0 . 002 , 1 ) 

9410 FORMAT(D1 6.8) 

CALL HEIGHT ( 0.14) 

CALL PHYSOR( XPHY , YP MV ) 

CALL AREA2D( AX , AY ) 

C 

C CLOSE GRAPH AND REOPEN FOR LABEL AXES NUMBER TYPES 

C 

CALL ENDGR( 0 ) 

CALL PHYSOR( XPHY , YPHY ) 

CALL AREA2D( AX , AY ) 

CALL XNONUM 
C 

C DETERMINE THE AXIS LABELS: 

C 

C X-AXIS LABEL AND NUMBER TYPE 

C 

CALL EXCMS( ’ CLRSCRN ’ ) 

WRITE( 4,9415) 

9415 FORMAT ( ’ DO YOU WANT INTEGERS ON THE X-AXIS? (Y/N )’,/) 

READ( 4 , 9400 ) ANS 
IF(ANS.EQ.YES) CALL XINTAX 

CALL XNAME( ’ A+NGLE! (?Q!, +DEGREES ! ) $ ’ ,100) 

C 

C Y-AXIS LABEL AND NUMBER TYPE 

C 

VRITE( 4,9417) 

9417 FORMAT ( ’ CHOOSE THE Y VARIABLE TO BE PLOTTED (II):',/, 

X ’ 1 V (WATT) ’ , / , 

X ’ 2 DJ ’ , / ) 

READ( 4 , * ) IYVAR 
IF (IYVAR.EQ.2) THEN 

CALL YNAME( ’D+J ! $ * ,100) 

ELSE 

CALL YNAME(’W ( W+ATT ! )$ ’ , 1 00 ) 

ENDIF 

C 

C ROTATE NUMBERS ON Y-AXIS AND REVERSE TICKS 
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c 

CALL YAXANG( 0 . ) 

CALL XREVTK 
CALL YREVTK 
WRITE( 4 , 7900 ) 

7900 FORMAT( ’ INPUT YMIN OF THE FIRST SUBPLOT (FI 6. 8 ):’,/) 
READ (4,*) YMIN 
CALL YAXEND (’ NOLAST ’ ) 

YMAX = YMIN+ . 1 

CALL GRAF( 0 . , 1 . , 1 . , YMIN , YMAX , . 1 ) 

CALL THKFRM( . 027 ) 

CALL FRAME 
CALL ENDGR( 0 ) 

C 

C CLOSED GRAF FOR AXIS LABELS, NOW REOPEN TO PLOT 

C 

CALL PHYSOR(XPHY.YPHY) 

CALL AREA2D( AX ,AY ) 

CALL RESET ( ’ YAXEND ’ ) 

CALL RESET ( ' XNONUM ’ ) 

CALL YNAME( BLNK ,2 ) 

CALL XNAME(BLNK,2 ) 

CALL ENDGR( 0 ) 

C 

C POSITION OR REPOSITION ORIGIN AND SCALE PLOT 

C 

CALL PHYSOR(XPHY.YPHY) 

CALL AREA2D( AX , AY ) 

CALL RESET (’YAXEND 1 ) 

CALL YAXEND( ’NOFIRST’ ) 

C 

C DEFINE COORDINATE SYSTEM 

C KEEP OR CHANGE CONSTANTS 

C 

CALL EXCMS( ' CLRSCRN ’ ) 

WRITE( 4 , 9423 ) 

9423 FORMAT ( ’ INPUT XMIN: (FI 6. 8)’) 

READ ( 4 , * ) XMIN 

WRITE( 4 , 9424 ) 

9424 FORMAT( 1 INPUT XSTEP : (FI 6. 8)’) 

READ (4,*) XSTEP 

WRITE( 4 , 9425 ) 

9425 FORMAT ( 1 INPUT XMAX : (FI 6.8)’) 

READ( 4 , * ) XMAX 

WRITE( 4 , 9427 ) 

9427 FORMAT ( ’ INPUT YMIN: (FI 6.8)’) 

READ( 4 , * ) YMIN 
WRITE( 4 , 9426 ) 

9426 FORMAT ( ’ INPUT YSTEP : (FI 6.8)’) 

READ (4,*) YSTEP 

WRITE( 4 , 9429 ) 

9429 FORMAT( ’ INPUT YMAX: (FI 6.8)’) 
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READ(4,*) YMAX 
CALL EXCMS( ’CLRSCRN’ ) 

VRITE( 4 , 9330 ) 

9330 FORMAT ( ’ HOW MANY TICKS PER X-AXIS STEP? (12)’) 
READ( 4 , * ) ITICKX 

CALL XTICKS( ITICKX) 

WRITE( 4 , 9331 ) 

9331 FORMAT ( ’ HOW MANY TICKS PER Y-AXIS STEP? (12)’) 
READ( 4 , * ) ITICKY 

CALL YTICKS( ITICKY) 

IF( XSTEP . EQ . 0 . 0 ) XSTEP = ( XMAX-XMIN ) / AX 
IF ( YSTEP . EQ . 0 . 0 ) YSTEP = ( YMAX - YMIN ) / AY 
CALL GRAF( XMIN , XSTEP , XMAX , YMIN , YSTEP , YMAX ) 

CALL YNONUM 

CALL YGRAXS( YMIN, YSTEP, YMAX, AY , ’ $ ’ , - 1 00 , AX , 0 . 0 ) 
CALL RESET (’YNONUM') 

CALL XNONUM 

CALL XGRAXS( XMIN, XSTEP, XMAX, AX, ’ $ > , - 1 00 , 0 . 0 , AY ) 
CALL RESET (’XNONUM 1 ) 

C 

C CURVE PLOTTING 

C 

CALL LINESP( 1 .8) 

CALL SPLINE 
CALL THKCRV (0.02) 

CALL GRACE (0.00) 

C 

C DRAW A CONNECTED 

C 

C MULTIPLE CURVE RETURN POINT 
2000 CONTINUE 

CALL EXCMS( ’CLRSCRN’ ) 

WRITE( 4 , 9445 ) 

9445 FORMAT ( ’ YOU HAVE A CHOICE 

+ ’ 1 

+ ' 2 

+ ’ 3 

+ ’ 4 

+ ’ 5 

+ ’ PLEASE INPUT LINE 

READ( 4 , * ) ICT 
IF( ICT . EQ . 2 ) CALL DOT 
IF( ICT . EQ . 3 ) CALL DASH 
IF( ICT . EQ . 4 ) CALL CHNDOT 
IF( ICT . EQ . 5 ) CALL CHNDSH 
CALL EXCMS( ’CLRSCRN’ ) 

C 

C DATA INPUT FROM EXTERNAL DATA FILE 

C 

WRITE( 4 , 9455 ) 

9455 FORMAT ( ’ INPUT THE FILE DATA SET NUMBER (12)’) 
READ( 4 , * ) INN 



CURVE 



OF LINE TYPE: ’ ./ , 
SOLID LINE’ , / , 

DOTTED LINE’ , / , 

DASHED LINE’ , / , 

CHAIN DOT LINE’ , / , 
CHAIN DASH LINE’ ,/ , 
TYPE NUMBER: ( II ) ’ , / ) 
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NUCOLS = 2 
IN=0 

DO 2400 11 = 1 ,5000 

READ( INN, 9460) ( YX( II , JJ ) , JJ= 1 .NUCOLS) 

9460 FORMAT ( 4X , F8 . 3 » 2X , F20 . 18) 

IF( ( YX( I I , 1 ) .LE. 0.0). AND . ( I I . GT . 2 ) ) GO TO 8000 
IN= IN+1 

2400 CONTINUE 
8000 CONTINUE 
ICOL = 2 
DO 8010 11 = 1 , IN 
X( II )=YX( 11,1 ) 

Y(II)=YX(II, ICOL) 

8010 CONTINUE 

REWIND INN 
C 

C TERMINATE DATA READING & PLOT CURVE 

C 

CALL CURVE( X , Y , IN , 0 ) 

C 

C RESET LINE TYPE 

C 

IF( ICT . EQ . 2 ) CALL RESET(’DOT’) 

IF( ICT . EQ . 3 ) CALL RESET( ’ DASH ’ ) 

IF( ICT . EQ . 4 ) CALL RESET( ’ CHNDOT ’ ) 

IF( ICT . EQ . 5 ) CALL RESET (’ CHNDSH ’ ) 

CALL EXCMS( ’ CLRSCRN ’ ) 

WRITE( 4 , 9475 ) 

9473 FORMAT ( ' DO YOU WISH TO TERMINATE THIS SUBPLOT?’) 
READ( 4 , 9400 ) ANS 
IF( ANS.NE. YES) GO TO 2000 
C 

C END OF PLOT 

C WRITE ADDITIONAL LABELS 

C 

CALL ENDGR( 0 ) 

CALL AREA2D( AX , AY ) 

CALL EXCMS( ’CLRSCRN’ ) 

WRITE( 4,9480) 

9480 FORMAT ( ’ DO YOU WANT ADDITIONAL LABELS ?’,/) 

READ( 4 , 9400 ) ANS 

IF( ANS . EQ . ’ N ’ ) GO TO 9900 

WRITE( 4 , 9483 ) 

9483 FORMAT( ’ YOU HAVE A CHOICE OF FUNCTION :’,/, 



+ > i GAUSSIAN FUNCTION’,/, 

+ • 2 LEVEL FUNCTION’,/, 

+ > 6 TRAPEZOIDAL FUNCTION’,/, 



+ ’ INPUT THE NUMBER OF FUNCTION TYPE (II) :’,/) 

READ (4,9401 ) NF 
IF (NF.EQ. 1 ) THEN 

INF = ’GAUSSIAN FUNCTIONS’ 

WRITE (4,9484) 



81 



9484 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 



FORMAT ( ’SPECIFY BUNCH LENGTH (F4.2 

> i 0.24 CM ’ , / , 

- 2 0.88 CM 

> 3 i.20 CM ’ 

> 4 i.76 CM ’ , / , 

> 5 2.40 CM ’ ,/, 

' 6 2.63 CM ’ ,/ , 

> 7 3.50 CM ’ , / , 

> 8 4.38 CM 

’INPUT THE NUMBER ( II ) : ’,/) 

READ( 4 , 9401 ) NBL 
IF (NBL.EQ. 1 ) THEN 

IBL= ’BL = 0.24 +CM ! $ ’ 

ELSEIF( NBL.EQ. 2) THEN 

IBL= ’BL = 0.88 +CM ! $ ’ 

ELSEIF( NBL.EQ. 3) THEN 

IBL= ’BL = 1.20 +CM ! $ ’ 

ELSEIF( NBL.EQ. 4) THEN 

IBL= ’BL = 1.76 +CM ! $ ’ 

ELSEIF( NBL.EQ. 5) THEN 

IBL= ’BL = 2.40 +CM! $ ’ 

ELSEIF( NBL.EQ. 6) THEN 

IBL= ’BL = 2.63 +CM!$’ 

ELSEIF( NBL.EQ. 7) THEN 

IBL= ’BL = 3-50 +CM ! $ ’ 

ELSE 

IBL= ’BL = 4.38 +CM!S’ 



9486 



ENDIF 

ELSEIF(NF.EQ.2) THEN 
INF= ’LEVEL FUNCTIONS’ 

VRITE( 4 , 9486 ) 

FORMAT( ’SPECIFY LEVEL FUNCTION LENGTH (F4.2 

+ ’ 1 1.00 CM ’ , / , 

+ ’ 2 3.00 CM ’ , / , 

+ ’ 3 4.50 CM ’ ,/, 

+ ’ 4 5.00 CM ’ ,/, 

+ ’ 5 6.00 CM ’ , / , 

+ ’ 6 7.00 CM ’ ,/, 

+ ’ 7 8.00 CM ’ ,/, 

+ ’ 8 9.00 CM ’ ,/, 

+ ’INPUT THE NUMBER( II ): ’,/) 

READ( 4 , 9401 ) LF 
IF (LF.EQ. 1 ) THEN 

LFL= ’LFL = 1.00 +CM! $ ’ 

ELSEIF( LF . EQ . 2 ) THEN 

LFL= ’LFL =3.00 +CM!$’ 

ELSEIF( LF . EQ . 3 ) THEN 

LFL= ’LFL =4.50 +CM!$’ 

ELSEIF( LF . EQ . 4 ) THEN 

LFL= ’LFL = 5.00 +CM!$’ 

ELSEIF(LF . EQ . 5 ) THEN 

LFL= ’LFL = 6.00 +CM!$’ 
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ELSEIF( LF . EQ . 6 ) THEN 

LFL= ’ LFL =7.00 +CM!$’ 
ELSEIF( LF . EQ . 7 ) THEN 

LFL= ’LFL = 8.00 +CM!$’ 
ELSE 

LFL= ’LFL = 9.00 +CM!$’ 



9501 



ENDIF 

ELSEIF( NF. EQ . 6 ) THEN 

INF= 'TRAPEZOIDAL FUNCTIONS 
WRITE( 4,9501 ) 

FORMAT ( ’SPECIFY 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 



9502 



INPUT THE 
READ( 4 , 9401 ) TLL 
IF (TLL.EQ. 1 ) THEN 
ITLL= ’ TFTL = 0.07 
ELSEIF (TLL.EQ. 2) 
ITLL= 'TFTL = 2 
ELSEIF (TLL.EQ. 5) 
ITLL= 'TFTL = 3 
ELSEIF( TLL.EQ. 4) 



TOP LEVEL LENGTH : ’ 

1 0.07 CM ’ ,/. 

2 2.94 CM ’ , / , 

3 3.185 CM ’ ,/ 

4 4.0 CM ’ , / , 

5 5.0 CM ’,/, 

6 5.04 CM ’ ,/, 

7 8.96 CM ’ ,/ , 

NUMBER (II): 



./ 



+CM ! $ ’ 
THEN 

94 +CM ! $ ’ 
THEN 

185 +CM ! S ' 
THEN 



ITLL= ’TFTL = 4.0 +CM!$' 
ELSEIF( TLL.EQ. 5) THEN 



ITLL= ’TFTL 
ELSEIF( TLL.EQ 
ITLL= ’TFTL 
ELSE 

ITLL= ’TFTL 
ENDIF 

WRITE( 4,9502) 
FORMAT ( ’SPECIFY 



= 5.0 +CM ! $ ’ 
6 ) THEN 
= 5.04 +CM ! $ ' 

= 8.96 +CM ! $ ' 



+ 

+ 

+ 

+ 

+ 



,/) 



BASE LENGTH : ’ , / , 

1 5.215 CM ’ ,/, 

2 7.0 CM ’,/, 

3 8.0 CM ’,/, 

4 8.96 CM ’ ,/, 

INPUT THE NUMBER (II): 

READ( 4 , 9401 ) BL 
IF (BL.EQ.1 ) THEN 

IBL= ’ TFBL = 5.215 +CM!S’ 

ELSEIF( BL . EQ . 2 ) THEN 

IBL= ’TFBL =7.0 +CM! S' 

ELSEIF( BL . EQ . 3 ) THEN 

IBL= ’TFBL = 8.0 +CM!S’ 

ELSE 

IBL= 'TFTL = 8.96 +CM’S’ 

ENDIF 



./) 
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ENDIF 

WRITE( 4 , 9485 ) 

9485 FORMAT ( ' SPECIFY RADIATION MODE (II):’,//, 

+ ’ 1 CERENKOV ’ , / , 

+ ’ 2 — - SUB-CERENKOV’,/, 

+ ’ INPUT THE NUMBER: ’,/) 

READ( 4 , 940 1 ) NMOD 
IF( NMOD . EQ . 1 ) THEN 

IMOD = ’ C+ERENKOV ! $ 1 

ELSE 



IMOD = ’ +SUB- ! C-t-ERENKOV ! $ 1 
ENDIF 

WRITE( 4 , 9487 ) 

9487 FORMAT ( ’ SPECIFY ENERGY OR NBETA (II): 
X ’ FOR E ’ , / , 

X ’ 1 — E = 2.555 MEV’ , / , 

X ’ FOR NBETA’ . / , 

X ’ 2 --- NBETA = 0.99 NM’ ,/ , 

X ’ 3 -- NBETA = 0.98 NM’ ,/ , 

X ’ 4 -- NBETA = 0.97 NM’ , / , 

X ’ 5 — NBETA = 0.96 NM’ ,/ , 

X ’ 6 -- NBETA = 0.95 NM’ ,/, 

X ’ INPUT THE NUMBER: ’,/) 

READ( 4 , 9401 ) NB 
IF( NB . EQ . 1 ) THEN 

C INB = ’ +N ! %B ! = 1 .00$’ 

INB = ’E = 2.555 M+E!V$’ 

ELSEIF( NB . EQ . 2 ) THEN 

INB = ’+N!*B! = 0.99$’ 

ELSEIF( NB . EQ . 3 ) THEN 

INB = ’ +N ! ! = 0.98$’ 

ELSEIF( NB . EQ . 4 ) THEN 

INB = ’ +N ! ! = 0.97$’ 

ELSEIF(NB . EQ . 5 ) THEN 

INB = ’+N!%B! = 0.96$’ 

ELSE 



INB = ' +N ! ! = 0.95$’ 

ENDIF 

CALL EXCMS( ’CLRSCRN’ ) 

WRITE( 4 , 9488 ) 

9488 FORMAT ( ’ SPECIFY ETA (II):’.//, 

X ’ 1 — ETA = 12 ’ , / , 

X ’ 2 -- ETA = 24 ’ , / , 

X ' 3 — ETA = 37 ’ , / , 

X ’ 4 -- ETA = 74 ’ , / , 

X ’ INPUT THE NUMBER: ’,/) 

READ( 4 , 9401 ) NETA 
IF( NETA . EQ . 1 ) THEN 
IETA = ’ %C ! = 12 $’ 

ELSEIF( NETA . EQ . 2 ) THEN 
IETA = ’ %C ! = 24 $’ 

ELSE IF ( NETA . EQ . 3 ) THEN 



,//, 
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9490 



IETA = 
ELSE 
IETA = 
END IF 
WRITE( 4 
FORMAT ( 



’%C< = 37 $’ 



' %C ! = 74 $’ 



9490) 

SPECIFY 



BEAM PATH LENGTH (II):’,//, 



1 -- L = 21 CM (ETA= 1 2 ,N=6 ) ’ , / , 

2 -- L = 42 CM (ETA=1 2 ,N=3 ) ’ , / , 

3 -- L = 64.8 CM ( ETA=37 , N= 6 ) 1 , / , 

4 -- L = 126.2 CM ( ETA= 1 2 , N= 1 ) 1 , / 

5 -- L = 129.5 CM (ETA = 37,N=3 ) ’ , / 



INPUT THE NUMBER 
READ( 4 , 9401 ) NL 
IF( NL . EQ . 1 ) THEN 

INL = 'L = 21 +CM ! $ ’ 
ELSEIF( NL . EQ . 2 ) THEN 
INL = ’L = 42 +CM ! $ ’ 
ELSEIF( NL . EQ . 3 ) THEN 
INL = ’L = 64.8 +CM ! $ ’ 
ELSEIF( NL .EQ . 4 ) THEN 



./) 



9492 



INL = 
ELSE 
INL = 
ENDIF 
WRITE(4 
FORMAT ( 



’L = 126.2 +CM ! $ 1 



’L = 129-5 +CM ! $ ' 



9492) 

SPECIFY 
X ’ 1 -- 

X ’ 2 -- 

X ’ 3 — 

X ’ INPUT THE 

READ( 4 , 9401 ) NRAM 
IF(NRAM.EQ. 1 ) THEN 
IRAM = '*L! = 10.52 



WAVE LENGTH 
RAMDA = 
RAMDA = 
RAMDA = 
NUMBER: 



(II ): * .//. 

10.52 CM (N=1 )’ ,/ 
3.50 CM ( N=3 ) ’ ,/, 
1 .75 CM ( N = 6 ) ’ ,/, 
’./) 



+CM! $ ' 



ELSEIF( NRAM. EQ . 2 ) THEN 



IRAM = 
ELSE 
IRAM = 
ENDIF 



’ %L ! 
’ %L ! 



=3.50 +CM ! $ ’ 



= 1.75 +CM ! $ ' 



XMES = 2.880 
YMES = 6.130 
IMESS = IMOD 

CALL MESSAG( IMESS, 100, XMES, YMES) 



YMES = YMES - 0.250 
IMESS = INB 

CALL MESSAG( IMESS, 100, XMES, YMES) 



YMES = YMES - 0.250 
IMESS = IETA 

CALL MESSAG( IMESS, 100, XMES, YMES) 
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YMES = YMES - 0.250 
IMESS = INL 

CALL MESSAG( IMESS, 100, XMES, YMES) 

C 

YMES = YMES - 0.250 
IMESS = IRAM 

CALL MESSAG( IMESS, 1 00, XMES, YMES) 

C 

YMES = YMES - 0.250 
IF( NF . EQ . 1 ) THEN 
IMESS= IBL 

ELSEIF( NF . EQ . 2 ) THEN 
IMESS=LFL 

ELSEIF( NF . EQ . 6 ) THEN 
IMESS= ITLL 

CALL MESSAG( IMESS, 100, XMES, YMES) 

YMES = YMES- 0.250 
IMESS = IBL 
ENDIF 

CALL MESSAG( IMESS, 100, XMES, YMES) 

C 

IF(HH.EQ.fi) THEN 
XMES = 1.20 
ELSE 

XMES =1.50 
ENDIF 

YMES = -2.25 
IMESS = INF 

CALL MESS AG( IMESS, 1 00, XMES, YMES) 

C 

9900 CONTINUE 

CALL EXCMS( * CLRSCRN ’ ) 

IF( IOUT . EQ . 2 ) THEN 
WRITE( 4 , 9493 ) 

9493 FORMAT ( 1 < CLEAR > , THEN READY FOR HARD COPY') 
ENDIF 

CALL ENDPL( 0 ) 

WRITE( 4 , 9495 ) 

9495 FORMAT ( * DO YOU WANT TO MAKE ANOTHER FIGURE?',/) 
READ( 4 , 9400 ) ANS 
IF( ANS.NE.YES) THEN 
GO TO 9999 
ELSE 
ENDIF 

CALL RESET ( 'ALL' ) 

GO TO 100 

9999 CONTINUE 

CALL DONEPL 

STOP 

END 
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